Threshold and high-frequency behavior of dipole-bound anion 
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An explicit analytic description is given for dipole-bound anion (DBA) as an excess electron bound 
to the molecular neutral due to its dipole moment. The calculated DBA photodetachment cross- 
section displays oc lo~ 2 behavior for large ui, in complete accordance with the experimental data 
[Bailey et ai, J. Chem. Phys 104, 6976 (f996)]. At the threshold the photodetachment cross-section 
displays the Gailitis-Damburg oscillations. 

PACS numbers: 33.80.Eh, 33.15.Ry 
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Considerable attention is payed presently to so- 
called dipole-bound anions (DBA), i. e. molecular 
negative ions in which the excess electron is bound 
to the molecular neutral (MN) due to its dipole mo- 
ment . Such structures play important role in var- 
ious physical, chemical and biological processes |2|. 

As early as in 1947 Fermi and Teller [3j in their 
analysis of meson capture by hydrogen atoms noted 
that a fixed point dipole d > 1.625 D can bind 
an electron to an infinitely large number of bound 
states. A number of subsequent studies taking into 
account the finite dipole effects presence of a 
short-range repulsive core potential |5j, rotational 
excitation , polarization and quadrupolar interac- 
tion (see Ref. |7j and references wherein) yield the 
minimum dipole moment 2-2.5 D required to form 
DBAs of common molecules 0, |U . 

DBAs for experimental studies are created mainly 
by free electron attachment under high pressure noz- 
zle expansion conditions 0,0,E3,O| and by charge 
transfer from Rydberg atoms [8J, [l3| . The created 
DBAs are investigated, e. o., by photoelectron spec- 
troscopy [a. ITo L lii T fia. I13I . 

Large-scale ab initio calculations (see, for in- 
stance, 113 E3- E3 an d references therein) of 
DBA structure studied the effects of correlation, or- 
bital relaxation, dispersion and charge-transfer in- 
teraction. However, despite the increasing accu- 
racy of such many-electron calculations, simplified 
one-electron local model potentials still have been 
proving their efficiency and suitability for large-scale 
computer simulations, as well as more analytical the- 
ories (^U and references wherein). 

In this Letter we use simple model which consid- 
ers DBA as an electron moving in a point dipole 
potential. The developed below simple analytic the- 
ory allows to explain some features of frequency- 
dependent DBA photodetachment. While the most 
of experimental photoelectron spectra are recorded 
at a constant photon frequency w, the Reference [14J 
reports the photodetachment rate oc w~ 2 . Such a be- 
havior is different from the photodetachment cross- 
section oc lo~'' 7 1" 1 in atoms and oc w~ 3 / 2 in atomic 
negative ions [2(J , and can be easily explained using 
the proposed analytical technique. Our model also 
predicts oscillatory behavior of the cross-section near 
the threshold. Such oscillations were predicted for 



electron scattered on hydrogen atom [21| and two- 
photon photodetachment of hydrogen anion with ex- 
citation of the residual atom [22T |. The oscillations 
(as well as the linear Stark effect) are due to constant 
dipole moment of excited states of nonrelativistic hy- 
drogen atom. Similar oscillations for electron scat- 
tering on polar molecules are discussed in [23j . 

We consider the MN to be a symmetric top whose 
rotational state is determined by its angular momen- 
tum j with the projection jq onto the molecule-fixed 
C-axis directed along the MN dipole moment d which 
is considered to be a point one. We also assume that 
MN remains in one of its vibrational states. The 
Hamiltonian of the excess electron moving in the 
field of the rotating MN is 

H = H Tot + T C + V(r, cos 0) = bj 2 + {b 6 - b c )j c 



2m r r 2 dr 



dr 



h 2 \ 2 P cos^ 



2m r 2 



(1) 



where 1 is the electron orbital momentum, b^ and b^ 
are the MN rotational constants, r is the electron 
radius-vector whose direction is determined by the 
spherical angles ip) in the molecule-fixed frame. 
Dimensionless dipole moment j3 — 2m c \e\d/h 2 = 
0.786 d(D), e and m e being the electron charge and 
mass. 

The solution of Schrodinger equation with the 
Hamiltonan can be found as a sum over the chan- 
nels corresponding to different MN states: 



(2) 



where the angular functions $j depend on -d, tp 
and on the MN coordinates. There are two limit- 
ing cases when the Schrodinger equation with the 
Hamiltonian is separable. The first case, Born- 
Oppenheimer approximation (BOA), takes place 
when the MN rotation is slow compared to the elec- 
tron motion. Quantitatively, the difference between 
the energy levels of the MN states with neighbor- 
ing j values in BOA is small compared with the 
difference between the excess electron levels. The 
opposite limiting case, inverse Born-Oppcnhcimcr 
approximation (IBOA) is fulfilled when the electron 
motion can be considered to be slow compared to the 
MN rotation. Thus in IBOA only one term remains 
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in j2J corresponding to the conserving MN angular 
moment j; such wavefunctions were used in rotation- 
ally adiabatic theory of rotational autodetachmcnt 
of DBA H3. 

The procedure of the separation of variables is 
given in detail in HUll (BOA) and [13 (IBOA). In 
both cases the expression © is effectively reduced 
to a single-channel radial function multiplied by the 
appropriate angular function. The latter has more 
simple form in BOA case, so all the following results 
(which are not sensitive to the angular dependence 
of the electron wavefunction) will be formulated for 
BOA case. 

The BOA angular functions Z = Z\ v satisfy the 
equation 



(-1 2 + /3 cos d)Zx v = r]Z Xr] 



(3) 



and can be decomposed over the familiar spherical 
harmonics: 

i 

Here the ^-projection A of the electron orbital mo- 
mentum is conserved due the axial symmetry. This 
non-spheric symmetry results in Z-mixing taken into 
account by the coefficients a^j . Together with these 
coefficients, Z\ v depend on the eigenvalue 77 of the 
operator 0, and r] — > 1(1+ 1), I = 0, 1, . . . for d — > 0. 
This three-diagonal eigenvalue problem was studied 
in [2^, |2(| for application to Rydberg states in po- 
lar molecules. The functions Z were used as early 
as in Debye's works on the Stark effect on polar 
molecules 28] . The same angular functions were 
used in analysis of the critical binding dipole mo- 
ment in DBA without studying the radial func- 
tions: for other references on usage of Z functions 
see [i^. 

The radial functions satisfy 



l__d 

r 2 dr 



dR 
dr 



-7T-R 



2m E 



n 2 



-i? = 0. 



(4) 



The solutions of the Eq. (0} for binding states (E = 
— h 2 y 2 /(2m c ) < 0) are McDonald functions (since 
they go to zero when r — > 00): 



R(r) 



--K p (xr), p = V»7 + 1/4- 



For real p these functions diverge at r — > 0: 



Kp(xr) ~ (ktY 



0. 



(5) 



so the existence (or absence) of binding electron 
states for 77 > —1/4 is determined by the behavior 
of the MN potential at small r. 

If we consider only the binding due to the MN 
dipole moment, then these states arise for 77 < —1/4 
when p is imaginary (it is easy to see that r\ = — 1/4 
corresponds to the above cited critical value d = 
1.625 D). Applying a formula similar to for imag- 
inary p we see that the functions R(r) oscillate for 



r — > according to JSJ) and R(r) ~ exp(— xr) when 
r — > 00. Assuming s = i\p\ the normalized radial 
wavefunctions are 



2sinh7rs\ 1 / 2 y. 



—=K is (yr), 
\ r 



- 1/4, r/<-l/4. 



(6) 



Oscillating behavior of wavefunctions make 
them different significantly from the wavefunctions 
of an electron in <5-potential. Interestingly, the func- 
tions © satisfy all the necessary boundary con- 
ditions for all E. In other words, a continuous 
spectrum of bound states arises in the field of a 
point dipole with 77 < —1/4. In this sense it is 
an unique case in quantum mechanics (excluding 
more trivial cases, e. g., quasicontinuous spectrum 
of bound states of an electron in a macroscopic po- 
tential box). The "falling to the center" (see the 
detailed comments in |29j) consists in that the ener- 
gies of the bound states are not limited from below 
and the wavefunctions corresponding to different en- 
ergies are not orthogonal. 

Of course, this continuous spectrum of bound 
states is interesting only from the mathematical 
point of view, and connected with the fact that the 
Sturm-Liouville problem J3J has a singularity at r = 
0. Physically, one should regularize the point dipole 
potential at small r. It can be achieved, for instance, 
by considering the MN as an extended dipole 
taking into account some short-range repulsive core 
potential 0, 0] . We use the simplest regularization 
model of the nonpenetrating core. The boundary 
condition in this model reads as R^ p (ro) = at 
some characteristic value of the MN radius fo- For 
yro <C 1 we can use the formula like (Jsj: 
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and obtain the electron energy spectrum: 



sin [s\n(yr /2) - argT(l + is)] = 0, 



E 



Xrjn 



2h 2 

5- exp 

m c r^ 



2ixn 1 „, 1 
+ -argr(l + is) }(7) 

s s 



where n = 1,2,... is the "principal" quantum num- 
ber and the dependence on 77 and A is included 
through s. Note that a formula similar to J2J de- 
scribes the spectrum of excited states of an elec- 
tron in the finite-size dipole field [3(| • The expres- 
sion J7J) should be applied carefully to low-excited 
DBA states; in the below calculation of photode- 
tachment rate we do not use this formula. 

For free electron states E — h 2 k 2 /(2m e ) > and 
the radial functions are expressed in terms of the 
Bessel functions whose order p may be cither real or 
imaginary: 

R k>s (r) = -i=[J ia (fcr) + J_ ia (fcr)], 77 < -1/4, 
2y/r 
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Rk, P {r) = -y=J p (kr), 77 > 



-1/4. 



(8) 



FIG. 1: Polar directional pattern of \Z\ V \ for A = and 
d = 3 D (77 = -0.74, dotted line), d = 4 D (77 = -1.18, 
dashed line), d — 5 D (77 = —1.66, solid line) 



Figure shows the polar directional patterns of 
\Z\ V \ 2 for cr-states with 77 < —1/4 at different dipole 
moment values. Qualitatively, these patterns resem- 
ble those obtained by ab initio calculations (see, 
e. g., 12]). The same patterns for 7r-states with 
77 > —1/4 are given in Fig. [3 

The cross-section of the photodetachment in 
dipole approximation is 
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FIG. 2: Polar directional pattern of \Z\ V \ 2 for A = 1 
and d = 3 D (77 = 2.36, dotted line), d = 4 D (77 = 2.5, 
dashed line), d = 5 D (77 = 2.6, solid line) 



Here a — e 2 /hc ~ 1/137, u> is the light field fre- 
quency, r s is projection of the electron coordinate 
onto the field polarization, ^ is the wavcfunction 

of the initial electron state. The wavefunction 
of the final state has the converging wave asymp- 
totics. Its angular part can be constructed using the 
technique developed in [25| for BOA and in [27J for 
IBOA. The frequency dependence of do~ e (x\i]; ui) is 
determined only by the squared radial integral Q 
between the functions 10 and (JSJ 



f {k/x)p' 2 F 1 {3/2 + p'/2 -is/2, 3/2 + p'/2 + is/2; 1 + p' --k 2 / x 2 ), 77' > -1/4 

^ K 1 (fc/^) is ' 2 F 1 (3/2 + i(s - s')/2, 3/2 + i(s + s')/2; 1 + s'; -fc 2 /^ 2 ) + (s — » s'), rf > -1/4 



(9) 



where k 2 / x 2 = hu/\E\-l and the Eq. 7.7(31) in ^ 
to rewrite the radial integrals in terms of Gaussian 
hypergeometric functions. Depending on the rela- 
tive energy of the MN rotation (i. e., BOA or IBOA), 
the parameters s and p' (s') depend on the dipole 
moment as well as on A, 77 and the MN rotational 
quantum numbers. The final p' (s') depends also on 
the radiation polarization e. 

The full cross section <j e {>c\r\\ ui) is the sum over 
the channels with different 77', similarly to the sum 
over the channels with different orbital momentum 
I in spherically symmetric case. The high-frequency 
limit k 2 / k 2 — > 00 can be obtained from using the 
asymptotics of the hypergeometric functions given 
in Eq. 2.1(17): 

a(ui) (x uj~ 2 , huj^\E\. (10) 

The dependence (|l()t is common for both perpen- 
dicular and parallel (e = x or e = z) polarization 
of the radiation as well as for all channels with 77 > 
— 1/4 and 77 < —1/4. This dependence differs from 
photodetachment cross section in s-states of atomic 
negative ions (o(ui) oc w -3 / 2 ) but agrees with the ex- 
perimental data [1JI for DBA. The difference from 



the zero-range potential model is due to the oc 1/ ^fr 
behavior of the wavefunctions JHJ in the small r do- 
main that is different form the behavior oc 1/r of 
the wavefunctions in the zero-range potential. The 
oc ui~ 2 behavior of <j{ui) holds also for the non- 
penetrating core model Q provided the de Broglic 
wavelength of the excess electron is larger than the 
effective core radius: k ~ ^J2m c uj/h < 1/tq. As- 
suming ro ~ 1 A, the photon frequencies ~ 3 eV 
used in Ref. [Ti| obey this condition. Actually, for 
much higher frequencies the ionization of the inner 
electronic shells becomes more probable than the de- 
tachment of the dipole-bound electron. 

From the Eq. one can also extract the behavior 
of the photodetachment cross-section at the thresh- 
old k 2 /x 2 -> 0. Note that the factor (k/x) p ' © 
suppresses the contribution of the rf > —1/4 chan- 
nels at the threshold, and the main contribution is 
given by the 77' < —1/4 channels. This fact has 
simple physical explanation. For 77' > —1/4, the 
centrifugal repulsion of the electron exceeds its at- 
traction by the dipole in the r ~ domain where 
the radial wavefunction is localized. Thus the ra- 
dial wavefunction is small in r ~ domain because 
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the electron with k — > cannot penetrate under the 
barrier. From the other hand, the barrier is absent 
for T}' < —1/4 so that the radial function does not 
disappear in r ~ domain and the contribution of 
the corresponding channels remains finite at k — > 0. 
So we obtain 



cos 



s ln(k/> 



■Ipa 



(11) 



where the phase ip as i does not depend on k. The re- 
sult 1)11(1 differs significantly from the power Wigner 
law of the threshold photoeffect cross-section but it 
is similar to the threshold behavior of reactions in- 
volving the excited hydrogen atom jEH |22| . Note 
also that the dependence (|ll|l cannot be obtained 
in Born approximation when the dipole influence on 
the free electron motion is neglected. 
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FIG. 3: Photodetachment cross-section in the near- 
threshold domain 

Figure shows the threshold behavior of the pho- 
todetachment cross-section u z when the radiation is 
polarized along the dipole direction. For the per- 
pendicular polarization, the cross-section tends to 
zero since for |A| > one has r\ > —1/4 (the first 
value J] < —1/4 arises for A = ±1 at d ~ 9.646 D). 
The cross-sections are given in relative units, their 
absolute values depend on the approximation (BOA 
or IBOA) used for the radial functions. It is seen 
that the cross-section is non-monotonic function of 
Hui/\E\ near the threshold. 
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